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ABSTRACT 


This  document  is  concerned  with  the  effect  of  atmospheric  turbulence  on  the 
tracking  error  of  a  millimeter  wave  radar  transmitting  over  a  propagation  path  near  the 
ground,  as  in  the  case  of  one  used  for  the  guidance  of  a  direct-fire,  surface-to-surface,  anti¬ 
tank  missile.  The  objective  is  to  estimate  the  likely  angle-of-arrival  error,  and,  based  on 
general  theory  and  specific  experimental  results  found  in  the  literature,  to  assess  the 
reliability  of  the  estimate. 

The  document  discusses  the  appropriateness  and  limitations  of  the  Rytov 
approximation  used  by  V.I.  Tatarsky  and  others  for  theoretical  calculations  relevant  to  the 
error  estimate.  It  also  discusses  the  applicability  of  certain  data  obtained  by  the  National 
Oceanic  and  Atmospheric  Administration  and  the  Georgia  Institute  of  Technology  in  a 
measurement  program  at  Flatville,  Illinois. 

The  primary  conclusion  is  that  the  angle-of-arrival  error  will  probably  be 
significant.  A  program  to  perform  measurements  in  a  greater  variety  of  environmental 
conditions  than  those  encountered  at  Flatville  and  to  include  direct  angle-of-arrival 
measurements  is  recommended. 
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I.  INTRODUCTION 


In  recent  years  the  possible  use  of  a  millimeter  wave  radar  for  guidance  of  a  direct- 
fire,  surface-to-surface,  anti-tank  missile  has  created  an  interest  in  the  effect  of  atmospheric 
turbulence  on  millimeter  wave  propagation  near  the  ground.  Despite  well  planned 
measurement  programs  (Refs.  1-2)  followed  by  analysis  of  the  resulting  data  (Refs.  3-4), 
based  on  much  experience  with  other  wavelengths  and  propagation  at  higher  altitudes 
(Refs.  5-7),  the  likely  contribution  of  turbulence  to  the  radar  tracking  error  remains 
uncertain. 

This  uncertainty  is  the  result  of  several  factors  associated  with  the  particular  nature 
of  the  tactical  requirements  imposed  on  the  system.  It  is  assumed  here  that  the  radar 
frequency  of  interest  is  94  GHz  and  that  the  propagation  path  is  parallel  to,  and  1  to 
2  meters  above,  the  ground.  It  is  also  assumed  that  target  ranges  of  interest  are  from  1  to 
5  km.  For  this  combination  of  parameters,  questions  arise  concerning  the  validity  of 
certain  approximations  used  (cf.  Ref.  5)  in  the  theory  of  propagation  in  turbulence.  In 
addition,  some  errors  and  misinterpretations  of  the  theory  in  prior  analyses  (Refs.  1-4) 
have  contributed  to  the  problem. 

The  objective  of  the  present  document  is  to  estimate  the  likely  angle-of-arrival  error 
based  on  data  reported  in  References  2  and  3  and  to  assess  the  reliability  of  the  estimate. 
The  conclusion  based  on  that  assessment  leads  to  a  recommendation  that  another 
experimental  program  be  conducted  under  a  greater  variety  of  temperature  and  humidity 
conditions  than  those  encountered  in  the  Flatville,  Dlinois,  experiment  (Ref.  2),  and  which 
would  include  direct  angle-of-arrival  measurements. 
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II.  SOME  DEFINITIONS  AND  BASIC  THEORY 


A.  LOCAL  HOMOGENEITY  AND  THE  STRUCTURE  FUNCTION 

Reference  5  and  its  updated  version  (Ref.  8)  by  the  same  author,  V.I.  Tatarski, 
provide  the  fundamental  theoretical  results  used  in  the  analyses  of  the  millimeter  wave 
turbulence  propagadon  experiments  mendoned  in  the  introducdon.  The  theory  is  based  on 
characterizing  the  turbulent  atmosphere  as  one  with  an  index  of  refracdon  that  is  a  spadal, 
i.e.,  three  dimensional,  random  process.  * 

In  some  applicadons  such  a  pnx:ess  may  be  stadsdcally  homogeneous,^  which,  if  it 
were  true  in  the  present  case,  would  imply  that  the  auto-correlation  of  the  index  of 
refracdon  is  a  funcdon  sadsfying 

n(^l)  “n(?2)  “Const.  ,  (1) 

where  the  overline  indicates  an  ensemble  average. 

If  the  process  is  also  isotropic  then 

K(f,-fj)  =  K(p)  , 

where 


However,  for  the  turbulent  atmosphere  the  index  of  refraction  is  not  homogeneous 
in  general.  It  is  more  often  the  case  that  the  index  is  locally  homogeneous,  which  is 
defined  by  a  weaker  condition,  valid  over  some  restricted  region. 


^  Reference  S  refers  to  the  spatial  random  process  as  a  random  field. 

^  The  term  "homogeneous"  is  a  generalization  to  three  dimensions  of  the  term  "stationary,"  which  refers 
to  a  one-dimensional  process. 
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[n(?i)-n(r,)]'=D„(7,-r2)  . 

The  first  equation  states  that  the  mean  of  the  difference  between  the  values  of  the  index  of 
refraction  at  two  points  is  a  function  of  the  displacement  vector  from  one  point  to  the  other. 
The  quantity  defined  on  the  left-hand  side  of  the  second  equation  is  called  a  structure 
function. 

In  particular, 

r 

|^n(L.yj,zJ-n(L,y2,Z2)J  =D„(p)  ,  (2) 

where 

That  is,  p  is  the  distance  between  two  points  lying  in  the  same  y,  z  plane,  defined  by  x  =  L. 
In  (2)  it  is  assumed  that  propagation  is  in  the  direction  of  increasing  x  and  that  the  structure 
function  is  defined  by  (2)  over  the  fixed  onhogonal  plane  at  x  =  L.  Reference  5  derives 
similar  structure  functions  for  the  phase  and  for  the  log  amplitude  fluctuations  of  a  plane 
wave  and  of  a  spherical  wave  propagating  through  a  turbulent  medium  characterized  by  an 
index  of  refraction  satisfying  a  relation  of  the  form  (2). 

The  turbulence  process  has  two  fundamental  parameters:  Iq,  which  is  the  inner 
scale  size  and  lower  bound  of  the  so-called  inertial  subrange  where  the  Kolmogoroff 
turbulence  model  predicts  a  dependence  (Ref.  6),  and  Lq,  which  is  the  outer  scale 
size  and  upper  bound  of  the  inenial  subrange.  Temperature  fluctuation  observations  near 
the  ground  suggest  that  the  quantity  Lq  is  roughly  one  third  the  height  of  the  observation 
point  above  ground  (Ref.  6).  Other  measurements  (Ref.  6)  indicate  that  the  quantity  Iq  is  of 
the  order  of  millimeters. 

B .  THE  WAVE  EQUATION  AND  THE  RYTOV  APPROXIMATION 
It  is  assumed  that  the  scalar  wave  equation 

V^u  +  k^n^  (?)  u  =  0  (3) 
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governs  the  electromagnetic  propagation,  where  n  is  the  index  of  refraction  for  the  assumed 
inhomogeneous  atmosphere  and,  because  the  inhomogeneities  are  small,  has  the  form 

n(r)  =  1  +  n^(r)  ,  I  n^C?)  I  «  1  . 

Tatarski  (Ref.  5,  Ch.  7)  uses  what  he  calls  "the  method  of  smooth  perturbations,"  but 
which  is  more  commonly  known  as  the  Rytov  method  (cf.  Ref.  7),  to  find  approximate 
solutions  of  a  pair  of  differential  equations  that  are  equivalent  to  (3), 

A^o  +  +  k^  =  0  ,  (4a) 

AVj  +  Vxjtj  •  ^2V\ir^  + +2k\j^rj^  =0  ,  (4b) 

where 

Vo  +  Vi=¥  > 

V  =  logu  , 

and  a  term  k^nj^  has  been  dropped  because  it  aJ  small  compared  to  the  term  2k2ni.  With 
the  assumption  that  yi  does  not  change  much  over  a  distance  of  the  order  of  a  wavelength, 
the  quadratic  term  in  the  gradient  of  y  i  can  also  be  dropped,  leaving 

Ay ^  +  2Vy^  •  Vy^  +  2k  n^Cr)  =  0 


which  provides  a  solution  to  equation  (4a)  to  use  in  (5),  and  neglecting  the  second 
derivative  of  yi  with  respect  to  x  in  (5),  leads  to  the  equation 

dy  2  ,  X 


The  differential  equation  (6)  embodies  the  Rytov  approximation,  sometimes  called  the 
parabolic  approximation  because  of  its  form. 

Reference  5  (p.  128)  shows  that  the  Rytov  approximation  is  valid,  given  the 
condition 
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L«lJ/X^  , 


(7) 


where  X  is  the  wavelength.  Obviously,  the  parameters  associated  with  the  systems  of 
interest  in  this  paper  will  not  satisfy  (7).  However,  it  is  clear  from  its  derivation  in 
Reference  5  that  the  condition  is  overly  restrictive.  The  derivation  of  the  Rytov 
approximation  starts  with  the  integral  representation 

.2  ikir-r'l 


dV'  , 


which  leads  to  the  approximate  relation  for  yi  that  satisfies  (6): 


/’‘Tm) 

f') - i - ^ - ^dV'  . 


where 


p=\/(y-y')  +(z-z')  . 


This  result  is  obtained  by  introducing  the  series  expansion 

ic|r-r1=k(x-x').^.-^.... 

and  dropping  all  terms  of  order  higher  than  quadratic  in  p.  Reference  5  obtains  the 
condition  (7)  by  treating  the  neglected  term  of  4th  order  as  an  estimate  of  the  error  and 
using  the  argument  that 

0 


X  -  x'  -  L  , 

in  the  important  region  of  integration.  In  deriving  (7)  Reference  5  does  not  take  into 
account  the  directionality  of  the  antenna  beam,  however,  a  narrow  beamwidth  would  justify 
dropping  terms  of  higher  order  in  powers  of  0L  without  the  necessity  of  imposing  the 


^  The  lateral  distance  p  is  the  chord  of  a  circular  arc  with  radius  L  and  subtending  an  angle  6. 


condition  (7).  An  argument  analogous  to  that  of  Reference  5  would  lead  to  the  condition 
with  the  antenna  aperture  diameter  D  replacing  Iq  in  (7): 


which  is  much  weaker  than  (7)  for  the  systems  of  interest  here. 

Another  justification  for  asserting  that  the  atmospheric  inhomogeneides  contributing 
to  the  forward  propagating  field  will  be  primarily  those  lying  in  a  narrow  cone  about  the 
forward  direction  is  the  general  principle  that  multiple  scattering  is  attenuated  in  a  random 
medium.  The  translator's  note  d  on  p.  265  of  Reference  5  gives  an  argument,  as  well  as  a 
reference  [I.  Kay  and  R.A.  Silverman,  "Multiple  Scattering  by  a  Stack  of  Dielectric  Slabs," 
Nuovo  Cimento,  V.  9,  Series  X,  Supplemento  No.  2,626  (1958)],  lending  support  to  this 
principle. 


In  terms  of  a  turbulent  atmospheric  index  of  refraction,  assumed  to  be  a  locally 
homogeneous  and  isotropic  random  variable.  Reference  5  uses  the  Rytov  approximation, 
i.e.,  a  solution  of  (6),  to  derive  the  log  amplitude  and  phase  structure  functions  for  a  wave 
propagating  through  the  medium.  From  these  results,  based  on  the  Kolmogorov-Obukhov 
model  of  turbulence  for  the  index  of  refraction  structure  function  (cf.  Ref.  5,  p.  31), 
Reference  5  derives  an  expression  for  the  variance  of  the  log  amplitude  %  of  the 
propagating  wave  (Ref.  5,  p.  153) 


where  is  the  structure  parameter  characterizing  the  index  of  refraction  structure  function 
in  the  Kolmogoroff  turbulence  model.  The  corresponding  log  intensity  variance  can  be 


obtained  from  ^8)  by  multiplying  the  right  hand  side  by  4. 


Experimental  evidence,  based  on  a  measurement  of  intensity,  indicates  (cf.  Ref.  9) 
that  the  Reference  5  prediction  of  the  log  amplitude  fluctuation  variance,  calculated  by 
means  of  the  Rytov  approximation,  saturates  (i.e.,  levels  off  and  ceases  to  grow)  with 
increasing  distance,  breaking  down  when,  according  to  Ref.  7,  the  Rytov  value  for  the  log 
intensity  variance  corresponding  to  (8), 

of  =  1.23k™  . 
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exceeds  0.3.^  On  the  other  hand.  Reference  7  shows  that  even  when 

(^  »  1 

the  structure  function  for  phase  fluctuations  changes  only  slightly. 


C .  THE  KOLMOGOROV-OBUKHOV  TURBULENCE  MODEL  AND 
ANGLE-OF-ARRIVAL  FLUCTUATIONS 


Reference  5  uses  the  Kolmogorov-Obukhov  model  of  turbulence  (cf.  Ref.  5,  p.  21) 
for  the  structure  function  of  the  index  of  refraction  (Ref.  5,  p.  150): 


D„(r)  = 


for  Ip  «  r  « 


for  r  «  Ip  , 


(9) 


which  implies  a  locally  isotropic  medium.  Associated  with  the  structure  function  Dn  is  the 
spectral  density  On;  they  are  related  by  (cf.  Ref.  5,  p.  47) 

PQ  OO  OO 

Dn(?)  =  2j  J  J  (l-cosi?  r)0„(i?)di?  . 

OO  OP  OO 

For  a  locally  isotropic  medium  this  becomes 

Dn(r)  =  87tj(l-^^)o„(K)K2dK  .  (10) 

0 

For  the  structure  function  given  by  (9)  Reference  5  suggests^  a  spectral  density 
given  by 

0^(K)  =  0.033C^k""^  ,  (11) 

where  it  is  assumed  that  On(K)  vanishes  for  values  of  k  larger  than  some  finite  value. 
Reference  8  (p.  76)  uses  a  more  precisely  defined  variation  of  (1 1): 


^  According  to  Reference  9,  the  threshold  for  saoiration  is  1.5  for  Oj  rather  than  0.3  for  However, 

Reference  9  states  that  saturation  occurs  for  propagation  distances  "as  little  as  250  m,"  while  Reference 
7  refers  to  saturation  for  prqiagation  distances  over  1.0  km. 

^  N.b.,  the  structure  function  given  by  (9)  is  defined  only  for  values  of  r  in  a  region  that,  as  is  evident 

from  the  restrictions  placed  on  r  in  (9),  is  not,  itself,  precisely  defined.  Therefore,  some  ambiguity 
exists  in  the  corresponding  spectral  density. 


where 


Reference  5  (pp.  136-137)  gives  the  equation 

DjCp)  =  4tcVl  J  1  -  Jq(kp)]  ^  1  +-^ sin  j  <D^  (K)ia/ic  (13) 

for  the  phase  structure  function  of  an  electromagnetic  wave  propagating  in  a  locally 
homogeneous,  isotropic  turbulent  medium.  In  addition.  Reference  5  (p.  226-227)  (cf.  also 
Refs.  6-8)  gives  the  equation 

- 2  Dib) 

(Aa)  =-^  (14) 

kb 


for  the  variance  of  the  corresponding  angle-of-arrival  fluctuations,  where  b  is  an 
interferometer  separation  or  a  telescope  aperture  diameter. 

Assuming  the  Obukhov-Kolmogorov  turbulence  law,  for  the  case  of  a  plane  wave. 
Reference  5  (p.  155)  gives 

D^(p)  =  2.91k^LC^p5^^  ,  (15) 


which  is  valid  when 


p  » , 


and  (p.  156) 

D^(p)=1.46k^LC^p5/3  , 


(16) 


which  is  valid  when 

1q  «  p  «  . 


For  the  circumstance  that  neither  of  these  conditions  is  satisfied.  Reference  5  (p.  155)  gives 
the  more  general  equation  (which  reduces  to  (15)  or  (16)  when  the  appropriate  condition  is 
satisfied) 


Dj(p)  =  2.91  LC^  -  0.62CJ  l‘ [  1  -  b^(p) ]  ,  (17) 

where 

bA(p)  =  BA(P)/BA(0) 

and 

Ba(P)  =  27t^(0.033)C^  ^’^P>  ( ^  “  • 

D.  PHASE  AND  ANGLE-OF- ARRIVAL  TEMPORAL  POWER  SPECTRUM 

According  to  the  so-called  frozen  turbulence  hypothesis  (attributed  by  Refs.  6  and  8 
to  G.  Taylor)  the  temporal  variadons  of  a  quantity  observed  at  a  point  in  the  atmosphere  can 
be  attributed  to  the  uniform  motion  of  the  atmosphere  past  the  point.  Assuming  the  frozen 
turbulence  hypothesis,  the  phase  structure  function  defined  by  (2)  becomes 

D^(p)  =  D^(v^x)  .  (18) 

where  vi  is  the  wind  velocity  component  transverse  to  the  electromagnetic  propagation 
direction  and  t  is  time.  According  to  Reference  6,  with  the  frozen  turbulence  assumption 
when  VI «  0,  it  is  justifiable  to  neglect  the  influence  of  random  motions  of  the  medium 
occurring  during  the  evolution  of  the  turbulenc*  only  if  the  condition 

Lq  » 

is  satisfied.^ 

Corresponding  to  the  temporal  function  defined  by  (18)  is  the  spectral  function 
W3(f)  defined  by 

oe 

Dj(v  X)  =  2  J  [  1  -  cos(2^x)]  Wj(f)  df  .  (19) 

0 

In  terms  of  the  quantity  defined  by 


^  However,  this  condition  will  not  usually  be  satisfied  in  the  cases  of  interest  here.  On  the  other  hand, 
according  to  Reference  10,  the  frozen  turbulence  model  may  still  be  valid  for  describing  the  high 
frequency  part  of  the  intensity  fluctuation  spectrum. 
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where 


f 

0  JlnXL  ’ 


Reference  8  (p.  268),  derives^ 


w^(f)  =  8.2- 10  ,  Cl«\  , 

Wj.(f)  =  4.M0“^C^k\v^^f“*^  .  n»l  .  (20) 

Using  (20),  Reference  8  also  derives  the  spectral  function  W5s(f)  of  the  spatial 
phase  difference  between  two  points  in  the  plane  x  =  L,  separated  by  the  distance  p.  The 
basic  relation  for  this  purpose  is  (p.  269) 


w  (f)  = 
5s 


from  which  Reference  8  obtains  (p.  271) 


=  4sin^(^) 


W,(f)  , 


Wg^(f)  =  0.033  cJk^Lv^^sin^^^^r*^  ,  p  » , 
w^^(f)  =  0.016  k^  sin^  j  f"*^  ,  p  «  /XL  , 

The  spectral  function  corresponding  to  the  angle-of-arrival  fluctuations  can  be 
obtained  from  W5s(f)  by  dividing  by  k2p2.  Figure  1  illustrates  the  curve  corresponding  to  a 
normalized  version  W(f)  of  the  spectral  function  defined  by  (22). 


2  The  numerical  coefficients  on  the  right-hariJ  sides  of  the  equations  in  (20)  differ  from  those  in 
Reference  6  by  a  factor  of  0.2S.  The  difference  a^iears  to  be  due  to  errors  introduced  into  the  Reference 
6  equations  (29)  and  (44),  in  which  a  factor  of  r/2  in  the  corresponding  equations  (IS)  and  (34)  of 
Reference  8  is  replaced  by  the  factor  2n. 
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III.  MEASUREMENTS  AND  DATA  ANALYSIS 


A.  THE  REDSTONE  ARSENAL  MEASUREMENTS 

Reference  1  describes  an  experiment  at  the  Redstone  Arsenal  airport  to  measure 
atmospheric  scintillation  effects  on  electromagnetic  wave  propagation  at  94  GHz  1.1m 
above  the  ground  over  a  range  of  1232  m.  To  measure  the  scintillation  at  optical 
wavelengths  the  authors  used  a  He-Neon  laser  over  a  range  of  300  m. 

The  authors,  who  apparently  assumed  that  the  structure  parameter  is  the  same  at 

2 

optical  as  at  millimeter  wave  frequencies,  observe  that  values  of  the  structure  parameter 
reported  in  the  literature  at  the  time  of  their  experiment  differed  widely,  ranging  from 
4  X  10“^^  m~^  to  7  X  10“*^  m“^^.  They  regard  the  higher  value  as  suspect  because  it  is 
so  much  higher  than  the  other  values.  However,  the  measurements  that  led  to  the  smallest 
value  were  made  50  m  above  the  ground,  whereas  those  leading  to  the  largest  were  made 
only  2  m  above  the  ground,  where  the  turbulence  is  normally  much  larger.® 

2 

As  for  their  own  estimate  of  C^;  the  optical  value  may  be  too  small  because  it  wai 

based  on  measurements  of  log  amplitude  fluctuations  over  a  distance  of  300  m  in  the  optical 

wavelength  region,  for  which  Reference  9  reports  that  saturation  can  be  expected  at  a 

distance  as  short  as  250  m.  The  effect  of  saturation  is  to  reduce  the  measured  intensity 

2 

fluctuations,  leading  to  an  underestimate  of  C^.  In  addition,  according  to  Reference  3,  in 
the  Flatville  experiment  the  structure  parameters  observed  at  the  optical  frequency  were  two 
orders  of  magnitude  smaller  than  those  observed  at  the  millimeter  wave  frequency. 

B .  THE  FLATVILLE,  ILLINOIS,  MEASUREMENTS 

Reference  2  describes  an  experiment  conducted  by  the  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  and  the  Georgia  Institute  of  Technology  (GIT)  near 
Flatville,  Illinois,  in  five  separate  data-taking  sessions:  (1)  June  3-July  15,  1983; 
(2)  October  30-December  11,  1983;  (3)  February  5-March  9,  1984;  (4)  May  29-July  5, 


®  Cf.  Reference  6. 
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1984;  (5)  January  29-March  26, 1985.  NOAA  collected  micrometeorological  data  and  GIT 
millimeter  wave  propagation  data. 

Propagation  observations  included  some  with  a  CO2  laser  beam  for  the  purpose  of 

2 

measuring  the  index  of  refraction  structure  parameter  as  well  as  the  wind  component 
orthogonal  to  the  propagation  path.  Both  the  optical  and  millimeter  wave  paths  were  at 
levels  about  3.7  meters  above  the  ground. 

The  NOAA  temperature  and  humidity  data  were  used  to  calculate  temperature  and 
humidity  structure  and  cross-structure  parameters,  which,  in  turn,  were  used  to  estimate 
for  both  the  optical  and  millimeter  wave  frequencies.  Reference  2  mentions  several 
problems  with  the  micrometeorological  measurements  during  the  first  session  as  well  as 
attempts  to  cure  the  problems  by  means  of  software. 

Reference  2  gives  an  example  of  log  intensity  fluctuation  variance  determined 

during  the  first  session,  in  July,  1983,  for  millimeter  wave  (173  GHz)  propagation  over  a 

2  2 

range  of  1.374  km  as  a  function  of  C^.  The  values  of  were  estimated  by  means  of  the 
micrometeorological  and  optical  data  at  six  different  times  of  the  day,  with  the  results,  as 
reported  in  Reference  2,  shown  in  Table  1. 


2 

Table  1.  Values  of  C  at  Various  Times  During  a  Day  In  July  1983  at 

Flatville,  Illinois  (Ref.  2) 


Time 

8:40  a.m. 

10:16  a.m. 

11:30  a.m. 

1 :40  p.m. 

3:50  p.m. 

6:50  p.m. 

B 

5.8 

5.9 

in 

4.7 

1.9 

Reference  2  presents  a  plot  of  the  log  intensity  fluctuation  variance^  obtained  from 

millimeter  wave  observations  at  the  same  times  versus  the  theoretical  values  of  the  variance 
2 

as  a  function  of  The  log  intensity  fluctuation  variance  calculations  are  based  on  the 
equation 


the  right  hand  side  of  which  is,  as  it  should  be,  four  times  that  given  in  Reference  5 
(p.  188)  f(x  the  log  amplitude  fluctuation  variance  of  a  spherical  wave. 


^  References  2  and  3  refer  incorrectly  to  the  quantity  plotted  as  intensity  instead  of  log  intensity. 
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However,  Reference  3,  discussing  the  same  experiment,  reports  somewhat 
2 

different  values  for  as  shown  in  Table  2.  Some  values  of  the  measured  intensity 
fluctuation  variance  also  appear  to  have  been  changed  in  the  Reference  3  plot  of  the 
measured  intensity  fluctuation  variance  versus  the  theoretical  based  on  C^j. 


2 

Table  2.  Values  of  C  at  Various  Times  During  a  Day  In  July  1983  at 

Flatville,  Illinois  (Ref.  3) 


Time 

8:40  a.m. 

10:16  a.m. 

1 1 :30  a.m. 

1:40  D.m. 

3:50  D.m. 

6:50  D.m. 

n 

5.3 

6.2 

5.2 

m 

1.9 

Since  the  authors  of  Reference  3  include  most  of  those  who  contributed  to 
Reference  2  and  the  publication  date  of  Reference  3  is  three  years  later  than  that  of 
Reference  2,  it  seems  reasonable  to  suppose  that  the  revisions  are  corrections  of 
presumptive  errors  in  the  original  work.  Perhaps  they  were  based  on  changes  in  the 
software  that  Reference  2  mentions  was  used  for  correcting  measurement  errors. 

The  abscissas  for  the  points  in  the  two  plots,  shown  in  Figs.  2  and  3,  are  the 

2 

theoretical  values  of  the  intensity  fluctuation  variance  based  on  the  values  of  shown  in 
Tables  1  and  2.  The  corresponding  Reference  2  ordinates  are:  0.(X)6,  0.013,  0.016, 
0.017, 0.018, 0.020.  For  Reference  3  they  arc  the  same  except  for  the  first  and  last,  which 
change  to  0.007  and  0.0225. 

A  linear  regression  Geast  square  error)  line  fit  to  the  points  in  the  Reference  2  plot 

has  a  slope  of  0.759  and  an  intercept  of  -3.83  x  10"^,  so  that  the  line,  in  effect,  passc«: 

through  the  origin  of  the  coordinate  system.  For  the  revised  Reference  3  points,  the 

regression  line  has  a  slope  of  0.796  and  an  intercept  of  -2.315  x  10"^,  so  that  the 

Reference  3  regression  line  passes  even  more  nearly  through  the  origin  than  the  Reference 

2  line.  Evidently,  in  either  case  the  measured  data  agree  quite  well  with  the  Reference  5 

2 

theory,  except  that  the  values  of  appear  to  have  been  underestimated  consistently. 

2 

The  explanation  for  the  low  (off  by  20  percent)  estimates  may  be  connected  with 
the  saturation  phenomenon  discussed  in  Reference  9:  experimental  data  indicate  that  for 
propagation  through  turbulence  at  some  optical  wavelengths  the  growth  of  the  log 
amplitude  fluctuation  standard  deviation  as  a  function  of  distance  levels  off  when 
1.5,  which  occurs  at  distances  as  small  as  250  m.  Reference  2  gives  1  km  as  the 
propagation  path  length  of  the  laser  beam  used  to  estimate  C^,  while  Reference  3  gives  the 
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revised  value  of  670  m  for  the  same  path  length.  In  either  case,  it  may  well  be  that  the 
propagation  distance  is  above  what  is  needed  for  saturation  to  occur,  even  allowing  for 
some  difference  in  the  CO2  laser  frequency  and  that  of  the  optical  beam  in  the  experiment 
discussed  in  Reference  9. 


0.5k^»  L”*  C* 

n 


Figure  2.  The  Normalized  Variance  of  Intensity,  Measured  at  Flatville, 
Plotted  as  a  Function  of  Its  Theoretical  Value.  The  solid  line 
shows  calculated  values.  The  dashed  line  Is 
a  least  square  fit  to  the  data. 
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Figure  3.  The  Normalized  Variance  of  Intensity  Versus  Its  lnertlal*Range 
Formula.  The  straight  line  shows  agreement  with  the  inertial  range  formula. 
The  symbol  gives  the  average  of  the  variances  from  antennas  1,  2,  and  3. 
The  vertical  separation  of  the  horizontal  bars  lying  above  and  below  each 
symbol  is  the  average  of  the  absolute  value  of  the  difference  between 
the  average  variance  and  the  variance  from  each  Individual  antenna. 
(Source  IEEE  Transactions  on  Geosclanca  and 
Ramote  Sensing, Mol  26,  No.  3,  May  1988) 
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from  one  fifth  to  one  third  the  height  of  the  propagation  path  above  the  ground,  while,  on 
the  other  hand,  according  to  Reference  2,  the  Flatville  millimeter  wave  path  was  only  about 
3.7  m  above  the  ground,  which  would  make  Lq  about  1.0  m. 

Reference  4,  using  the  values  =  5.9  x  10"^^  m“^^,  r  =  1  m,  and  a  hypothetical 

range  L  =  5000  m  in  (23),  finds  that  the  associated  standard  deviation  of  the  angle-of- 

arrival  fluctuations  is  126  prad,  and  that  extrapolating  the  Flatville  data  to  the  same 

conditions  gives  standard  deviations  ranging  from  85  to  200  |irad.  However,  assuming  the 

Flatville  frequency  of  173  GHz,  which  is  equivalent  to  a  wavelength  of  1.73  mm,  using  the 

same  parameter  values  in  (17),  and  calculating  the  angle-of-arrival  fluctuation  variance  by 

means  of  (14)  gives  a  different  range  of  angle-of-arrival  fluctuation  standard  deviations 

corresponding  to  the  Flatville  measurements.  This  is  evident  from  an  examination  of 

Table  3,  which  displays,  for  both  the  Flatville  range  of  1374  m  and  the  5000  m  range  used 

2 

in  the  Reference  4  example,  the  standard  deviations  calculated  this  way  using  the  values 
given  in  Table  2. 


Table  3.  Angle-of-^rrlval  Fluctuation  Standard  Deviations  Based  on  Table  2 

Values  for  p  s  1  m,  Vlfavelength  s  3  mm 


Time 

8:40  a.m. 

10:16  a.m. 

11:30  a.m. 

1:40  p.m. 

3:50  p.m. 

6:50  p.m. 

Angle  of  Arrival 

U1374  m 

125 

139 

150 

137 

123 

82 

oprad 

L=5000  m 

229 

254 

275 

252 

226 

152 

Reference  4,  having  remarked  that  the  fluctuations  observed  at  night  were  at  least  an 
order  of  magnitude  smaller  than  those  observed  during  the  day,  points  out  that  for 
millimeter  waves  the  worst  case  conditions  for  large  angle-of-arrival  fluctuations  should 
occur  on  a  hot,  humid  day.  An  example  cited  is  measurements  made  at  the  White  Sands 
Missile  Range  in  the  fall  of  1987,  a  few  days  after  the  occurrence  of  unseasonable  heavy 
rain.  Apparently,  larger  than  expected  fluctuations  of  94  and  140  GHz  waves  occurred. 
This  was  attributed  to  additional  contribution  to  the  turbulence  due  to  an  increased  humidity 
structure  constant  resulting  from  strong  sunlight  evaporating  the  water  in  the  wet  ground. 

Reference  4  concludes  that  amplitude  fluctuations  due  to  turbulence  should  not  be  a 
serious  problem  for  the  kind  of  millimeter  wave  radar  application  of  interest  but  that,  under 
worse  case  conditions,  angle-of-arrival  fluctuations  may  be.  Otherwise,  the  conclusion  is 
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that  in  most  circumstances  even  angle-of-anival  flucmations  should  be  less  than  a  serious 
problem. 

Reference  4  discusses  the  possibility  of  reducing  the  angle-of-arrival  fluctuations  by 
time  averaging.  However,  the  discussion  is  based  on  the  assumption  that  the  statistical 
behavior  of  the  angle-of-arrival  depends  on  phase  fluctuation  statistics,  which,  as  remarked 
earlier,  is  erroneous. 
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IV.  CONCLUSIONS 


This  document  argues  that  Reference  2  underestimates  certain  index  of  refraction 
structure  parameters  derived  from  the  Flatville  measurements  and  that  the  Reference  4 
method  of  calculating  angle-of-arrival  fluctuation  statistics  is  incorrect  theoretically. 
Nevertheless,  the  conclusion  of  Reference  4  that  millimetei  wave  amplitude  fluctuations 
should  have  a  negligible  effect  on  the  performance  of  a  direct  fire,  anti-tank  missile  system 
appears  to  be  correct.  On  the  other  hand,  the  conclusion  of  this  document  concerning 
angle-of-arrival  fluctuations  is  somewhat  more  pessimistic  than  that  of  Reference  4,  which 
predicts  that  they  will  have  an  important  effect  on  tracking  accuracy  only  in  hot,  humid 
weatlier.  Tliis  conclusion  is  based  on  its  estimate  of  angle-of-arrival  values  corresponding 
to  the  Ratville  structure  parameter  measurements  and  the  observation  that  those  values  are 
large  only  near  midday  when  both  the  temperature  and  humidity  were  high.  However,  the 
angle-of-arrival  standard  deviation  values  given  in  Table  3  of  the  present  document,  at  least 
for  the  5000  m  tange,  are  unacceptably  large  throughout  the  entire  day. 

Reference  4  recommends  a  measurement  program  to  get  values  of  the  index  of 
refraction  structure  parameter  in  other  geographical  locations,  under  different 
environmental  conditions.  This  could  settle  the  question  of  how  much  atmospheric 
turbulence  will  affect  millimeter  wave  angle-of-arrival  fluctuations  in  other  than  hot,  humid 
weather,  providing  very  useful  data  pertinent  to  the  design  and  application  of  millimeter 
wave  radar  systems  having  propagation  paths  near  the  ground. 

But,  because  the  nature  of  the  index  of  refraction  spectrum  is  uncertain  at  the  low 
frequency  end  when  (cf.  Ref.  5,  p.  150) 

p  =  . 

and  because  of  the  possibility  of  saturation  at  some  ranges  of  interest,  perhaps  an  even 
more  useful  addition  to  the  program  would  be  to  accompany  the  measurements  of  intensity 
fluctuations,  like  those  done  at  Flatville  for  the  purpose  of  estimating  angle-of-arrival 
statistics  indirectly,  with  direct  angle-of-arrival  measurements.  Indeed,  Reference  10 
reports  on  just  such  an  experiment  involving  millimeter  waves  and  outlines  the  method 
used  for  measuring  angle-of-arrival  fluctuations. 
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